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Abstract: A library rich in CA/GT microsatellites was constructed from the Paralichthys olivaceus genome by com- 
bining biotin capture method and radioactive labeling hybridization. Five hundred and twenty six positive clones were ob- 
tained through twice screens. Sequencing confirmed 133 microsatellite loci (number of repeats 5) in 119 positive clones. 
Of these microsatellites, two (1.5%) had compound repeat motifs, 63 (47.37%) had perfect motifs and 68 (51.13%) 
had imperfect motifs. Primer pairs were designed in the flanking regions of 22 microsatelites and subjected to PCR amplifi- 
cation. In 8 artificial gynogenesis families, four pairs failed to amplification, one pair was monomorphic, and the rest were 
polymorphic with an average of 5.2 alleles per locus. Heterozygosities ranged between 0.375 and 0.846, PIC ranged be- 
tween 0.305 and 0.823. The results suggested that most of the microsatellites we isolated were qualified to be applied to 
the population genetic studies of P. olivaceus . 
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FS CA/GT WML Strid hh wie 


BH? IAAL, SAK, ARES ARR, mee 
a. PAKAR ERR BE BK OAT, WEL PARR 150070; 2. PUK ERE, OH OM 510275; 
3. LY KOR EMERE LE 200090; 4. Kék SR EES SRASER, WT KE 116023) 


HE: KAA WREAK SRN EA RAE AA OAR, TEM FF (Paralichthys olivaceus ) 
RAAPA AMAA CA/CT BERNAL SIPS. Pau, SEARS 526 TATE. WP PHY 119 
SE, FRR 133 SAAR ERE. RT TRA DBS (1.5%), SERB 63 4 (47.37%), FE 
TÆ 68 (51.13%). HITFFA RM 22 WHEE, M8 RA LHR RA RAD RBA BOOT. PCR 
SERRA, 4 S| WIT BRAT BHR E AMA, 1 ARRAES, ER 17 IAEE St, FH 
FES POE 5.2 NESE, RAEN 0.375 ~ 0.846, SAME BAHAY 0.305 ~ 0.823. RREH, ATEN 
KBD TB tic BE we AO BE AEE. 


Kid: FH; MEERE; WRAP; RRRA 
FARES: 0346.5 XRRR: A XES: 0254 - 5853(2005)06 - 0652 - 05 


Microsatellites, also known as simple sequence re- repetitive sequence motifs (Degnan & Arévalo, 2004). 
peats (SSRs), are regions of DNA that exhibit short Those motifs are often composed of 1 - 6 bp repeat se- 
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quences, such as CA, AGA, ATA and the like. In re- 
cent decades, microsatellites were widely employed in 
population genetic studies of numerous species, and 
such application is continuously expanding (Kohlmann 
et al, 2005; Mia et al, 2005). This drives the devel- 
opment of methods to isolate microsatellite from the 
genome. Microsatellites are featured cross species am- 
plification due to the conservativeness among relative 
species in their flanking region where primers can be 
designed (add the reference), which reduces the re- 
quirements of direct isolation in target species (Mia et 
al, 2005). However, direct isolation from genome, 
which avoid drawbacks of cross species amplification, 
such as allele dropout, null allele etc, is still the main 
source of microsatellites (Xu et al, 2001). Presently, 
there are mainly two ways to isolate microsatellites from 
genome directly, non-enriched and enriched method. 
The former is not only time-consuming (generally three 
to eight months) but also low net yield of microsatell- 
lites, while the latter is often time-saving (about one 
month) and high net yield of microsatellites (Dimsoski 
& Toth, 2001; Chang et al, 2005). Therefore, en- 
riched method gradually plays major role instead of 
non-enriched method. 

: Paralichthys olivaceus is one of the most impor- 
tantly cultured fish species around the oceans in North- 
ern China. Because the female grows faster than the 
male, artificial gynogenetic technology and fry femi- 
nization are regarded as effective strategies to approach 
monosexual reproduction and high productivity (Liu et 
al, unpublished). However, there are very limited re- 
ports concerning the molecular genetic analysis of artifi- 
cial gynogenesis in P. olivaceus . We isolated and char- 
acterized a set of microsatellites directly from the P. o- 
livaceus genome. 


1 Materials and Methods 


1.1 Microsatellite isolation 

A modified version of the protocol outlined by Car- 
leton et al (2002) was used in the construction of the 
microsatellite library of the P. olivaceus . 

Total genomic DNA from a single individual of P. 
olivaceus was digested with the restriction enzyme Sau 
3A I. Fragments ranging from 400 - 800 bp were ex- 
cised from a 1% low melting point agarose TAE gel and 
purified using QIAquick Gel Extraction Kit (Qiagen, 
CA). The double-stranded adaptor molecular A/B was 
prepared by mixing equal mole amount of oligonu- 
cleotides A (5'-GATCGTCGACGGTACCGAATTCT-3’) 
and B ( 5’-GTCAAGAATTCGGTACCGTCGAC-3’) , 


heating to 95 °C at 10 min, then slowly cooling down to 
10 © over a period of 4 h. The recovered DNA frag- 
ments were suspended in a 20 „L ligation reaction con- 
taining 6.6 ug of A/B double-stranded adaptors and 6 
Weiss units of T4 DNA ligase. Ligation was carried out 
at 16 C overnight. 

Excess adaptors were removed by washing on an 
Ultrafree column (Pall, USA) 2-3 times. Two micro- 
liter of the recovered adaptor-ligated fragments was used 
as templates to perform PCR in a volume of 20 pL 
(oligonucleotide B as primer) following the program: 
94 © 3 min; 10 cycles for 94 C 1 min, 58 C 1 min, 
72 CT 1 min; 72 C 10 min for final extension. The 
PCR products were also purified 3 times to discard the 
remaining primers, dNTPs and concentrated to an ap- 
propriate volume with Ultrafree column. 

For enrichment, the adaptor-ligated DNA frag- 
ments was denatured at 95 %C for 5 min, then hy- 
bridized to a biotinylated probe (5-Biotin-ATAGAATAT 
[CA ];s) in 50 uL hybridization solution (6 x SSC, 10 
pmol/L oligo B, 0.1% SDS, 0.3 pmol/L probe) at 68 
© for one hour. The DNA hybridized to the probe was 
separated and captured by streptavidin magnetic beads 
(Dynal Biotech ASA, Norway) at room temperature for 
20 min, followed by 4 washing steps, including twice 
in 6 x SSC/0.1% SDS at room temperature for 10 min, 
twice in 3 x SSC/0.1% SDS at 68 for 15 min, twice 
in 6 x SSC at room temperature, and twice in 0.1 x TE 
at room temperature quickly. The final step of 95 C for 
10 min in 50 pL 0.1 x TE is to denature the single- 
stranded DNA containing the target microsatellites from 
the beads. The purified single-stranded DNA was sub- 
jected to a second round of PCR according to the same 
procedure as the first round of PCR. PCR products 
were also recovered as mentioned above. 

The recovered PCR products were ligated into 
pGEM-T vector (Promega,USA) according to the man- 
ufacturer’s instructions. The ligation solution was used 
to transform the competent E. coli cells (DH5a), 
which then plated and cultured on LB agar plate con- 
taining 100 pg/mL ampicillin. The recombinants were 
picked out orderly to nitrocellulose membrane 
(Promega, USA) in order to identify positive clones by 
5'-[y-*P ATP labeled probe (CA) prior to sequenc- 
ing. 

1.2 Sequence analysis and primer design 

Positive clones were sequenced at Chinese Nation- 
al Human Genome Center (http://www. chgb. org. 
cn). Software VecScreen (http://www. ncbi. nlm. nih. 
gov/VecScreen/VecScreen. html) was used to remove 
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the sequences of T vector, and a computer program de- 
veloped by us was used to remove the adaptors. The se- 
quences containing motifs repeating more than 5 times 
were regarded as microsatellites. Twenty two pairs of 
PCR primers were designed in the flanking regions of 
the repeating motifs with the software package of Pre- 
mier Primer 5.0 and subjected to PCR amplification as 
follow.. 
1.3 Microsatellite genotyping of P. olivaceus by 
PCR 

The genomic DNA was obtained from eight moth- 
ers of artificially gynogenetic family of P. olivaceus . 
Amplification reactions were carried out in thermocycler 


of Taq polymerase. Thermal cycling conditions were as 
follows: 94 © 3 min; 40 cycles for 93 C 30 s, 46- 
59 C 30s, 72 © 30s; 72 C 10 min for final exten- 
sion (Tab. 1). A horizontal agarose gel (2%) and a 
vertical polyacrylamide gel (8%) electrophoresis sys- 
tem were used to genotype each of these microsatellites 
as reported (Sambrook et al, 2001; Xu et al, 2002). 
1.4 Data analysis 

The expected heterozygosity (He) and the poly- 
morphic information content ( PIC ) were calculated 
through a computer program written by us according to 
the following formulae (Nei, 1978; Liang et al, 
2004) . 



























































9 700 (Pekin Elmer, USA) in a 25 uL volume con- L2 
tainng 1 - 100 ng of template, 10 mmol/L Tris-HCl He =1- > Pi E1 
(pH 8.3), 50 mmol/L KCl, 200 pmol/L each dNTP, n=l al ai 
0.1% Triton X-100, 0.1% NP-40, 0.01% gelatin, PIC = 12) Pi- >) Di2P%P? (2) 
0.4 pmol/L of forward and reverse primers, and 1 unit E _ 
Tab. 1 Microsatellite primer pairs for Paralichthys olivaceus , the number of alleles, the heterozygosity, and the 
polymorphic information content 
Locus Motifs Sequence (5'-3') ee a He HG eee. 
Pac2HLJ CAi7 eee 52 5 0.753 1 0.7120 DQ097685 
Pac6HLJ CA34 eee ica AT 48 4 0.745 6 0.698 3 DQ097689 
Pacl4HLJ CA n A A 52 3 0.709 1 0.655 0 DQ097688 
Pac18HIJ CA6N,CA); PE ee mate ae 55 4 0.686 4 0.6374 DQ097687 
Pacl9HLJ CAx2N2CA3 a eae. 51 7 0.826 4 0.804 3 DQ097686 
Pac28HLJ TG23N2TG, n ee DE eS 46 2 0.50 0.375 DQ097690 
Pac34HLJ a ca, ee Š 59 8 0.845 7 0.826 8 DQ097694 
Pac38HLJ GT 3 i eee 51 5 0.753 1 0.7120 DQ097693 
F:CACATC CGCA 
Pac39HLJ CAg2 R aire ag wins A 52 4 0.704 1 0.649 9 DQ097692 
Pac44HLJ GTis PEE eet a 51 T 0.828 4 0.805 9 DQ097691 
CA28N2CA3 F:TT TGACA Tı 
PacSOHLJ Nols i R: ue 46 5 0.763 9 0.726 0 DQ097695 
PacSSHLJ CA2 te ee ec 51 5 0.764 4 0.727 3 DQ097699 
F: GA 
= FFTEGCACIOGACETECT Tg 
Pac88HLJ GTs2 EER ERARE RL 51 4 0.613 3 0.536 9 DQ097697 
Pac99HLJ TG22N2TGs EE ph ana 51 6 0.789 1 0.758 8 DQ097696 
F : CGTAGGCAGCTCCAGAGT 
Pac106HLJ GTi R: CCACAGOGACCTCAMACA 47 2 0.375 0.304 7 DQ097700 
F:TTCATTCT CT 
Pacl 1OHLJ CA\3N2CA3 R: PE RE 51 2 0.375 0.304 7 DQ097702 
Pac130HLJ TG. Ne6TGay eee 49 6 0.786 4 0.7514 DQ097701 
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Where He is the expected heterozygosity in a popula- 
tion; n is the number of alleles at one locus; p;, p; are 
the frequencies of ith, jth allele at one locus, j = i + 
1. 


2 Results 


2.1 Microsatellite isolation and characteristics 
Prior to double check with the radioactive labeling 
method, a total of 2000 clones in the first round screen- 
ing were obtained. The positive clones were identified 
by PCR method with the universal primers. Approxi- 
mately 100% clones were positive showed on the PCR 
results (Fig. 1). However, only 526 positive clones 
(26.3%) were found after the second round screening 
with the radioactive labeling method. Sequencing the 


119 positive clones confirmed that 111 positive clones 
contain 133 microsatellite loci (numbers of repeats > 
5), 79 of which have unique flanking sequences where 
primers were designed. Of these sequences, only two 
contained compound repeat motifs (1.5%); of the re- 
maining sequences there involved 63 perfect (47.37% ) 
and 68 imperfect (51.13%) repeat motifs. The most 
frequent motifs were CA/GT repeats (93.23%) as an- 
ticipated. Though the motif CT/GA was not screened, 
there was a small proportion (6.77% ) found in associ- 
ation with the target motifs. 

Fig. 2 showed that the repeat length of CA/GT mi- 
crosatellite loci frequently occurred in 11 - 30, and 
when the length of repeat motif increased, the number 
of microsatellite loci was on the decline. 
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Fig. 2 Numbers of CA/GT microsatellite loci in 
accordance with repeat length 


2.2 Genetic analysis of 18 microsatellite loci in 
P. olivaceus 
Among 22 primer pairs designed at microsatellite 
loci, 18 of which amplified the target bands from the 
genomic DNA of 8 female P. Locus 
Pac78HLJ was proved monomorphic and all other loci 


olivaceus . 


were polymorphic. The data analysis showed that the 
number of alleles per locus averaged 5.2, the heterozy- 
gosity of each locus ranged between 0.375 and 0.846, 
and PIC ranged between 0.305 and 0.823 (Tab. 1). 


3 Discussions 


A variety of methods to isolate microsatellites were 
reported. However, the enriched method was the most 


Ten positive clones obtained in the first round screening 
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frequently employed by many research groups (Zhang 
et al, 2002; Chang et al, 2005). In our lab, mi- 
crosatellites have been isolated from more than 20 
aquaculture species by the biotin-capture method, and 
detailed information of these studies is available at our 
web site, www. fishbreeding. org. 

Theoretically, all of the clones contain the target 
microsatellites by the method of biotin capture (Reddy 
et al, 2001). However, only 10% - 90% of the 
clones are positive according to the results obtained in 
our lab and some other groups (Zhang et al, 2002; 
Chang et al, 2005). As shown in Fig. 1, if simply 
screening the library with the biotin capture method, 
the positive rate of clones was 100%, whereas only 
26.3% were really positive through the double check 
with radioactive-labeled probe. This may be due to 
three reasons: Firstly, after DNA hybridized to bi- 
otinlated probe, a series of subsequent washing 
steps didn’t remove the non-specific fragments com- 
pletely. Secondly, the streptavidin-coated magnetic 
beads have been stored too long (more than six 
months) before they were used, which seemed a key 
factor impacting the positive rate (Reddy et al, 2001). 
Finally, the screen with the radioactive-labeled probe 
led to the loss of some small fragments containing mi- 
crosatellites, and only a single set of oligos restricted 
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the potential range of microsatellite species being de- 
tected (Chang et al, 2005; Degnan & Arévalo, 
2004) . 

Of 22 microsatellite loci, 17 were highly polymor- 
phic, with an average of 5.2 alleles per locus, het- 
erozygosities ranging between 0.375 and 0. 846, and 
PIC varying between 0.305 and 0.823. The results 
indicated that a majority of the loci isolated in our re- 
search were potentially valuable for population genetic 
studies, marker densification on genetic linkage map 
(Maria et al, 2003), and marker-assisted selection 
(MAS) in gynogenetic families. 

As shown in Fig. 2, the vast majority of mi- 
crosatellite alleles were less 45 repeat units in length, 
which was similar to the results reported (Carleton et 
al, 2002; Webster et al, 2002; Garza et al, 1995). 
This may be caused during the evolutionary course of 
microsatellite by factors such as selection, biased mu- 
tation, gene conversion, and random mutation in re- 
peat units (Garza et al, 1995; Zhang & Zhang, 
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